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ABSTRACT 

Using a Monte Carlo technique, we investigate the center shifts that are expected 
to occur for the broad Ha wings and the He II A 6545 feature as a function of the 
neutral hydrogen column density in symbiotic stars. These two features are proposed 
to be formed via Raman scattering of UV continuum around Ly/3 and He II A 1025 
emission line by neutral hydrogen. The strengths of these two features are determined 
by a combination of various physical parameters including the covering factor and the 
neutral hydrogen column density (Nhi) of the scattering region. The branching ratio of 
Raman scattering to Rayleigh scattering for UV radiation around Ly/3 is a non-linearly 
increasing function of the wavelength, which results in enhanced Raman optical fluxes 
redward of the Ha line center as Nhi increases. However, we find that the amount 
of wing center shift is quite small about 20 km s _1 as Nhi increases from 10 20 cm~ 2 
to 10 21 cm -2 . Assuming that He II A 1025 emission is characterized by a Gaussian 
profile, the Raman scattered He II 6545 feature exhibits near Gaussian profiles with 
the peak shifted redward for Nhi < 10 22 cm~ 2 . The redward center shift amounts to 
1 A for Nhi ~ 10 20 cm~ 2 and decreases as Nhi increases up to Nhi ~ 10 22 cm~ 2 , 
above which no center shift is observed. The redward peak shift is due to the fact 
that the incident emission profile is symmetric with respect to He II A 1025 line center 
whereas the Raman conversion rate is increasing towards the line center of Ly/3. We 
emphasize that the determination of Nhi by locating the exact peak position of the 
He II 6545 feature will lift the degeneracy to allow a more accurate estimate of the 
covering factor of the neutral region, providing strong constraints on the mass loss 
process occurring in symbiotic stars. 

Key words: scattering — radiative transfer — binaries: symbiotic — mass-loss — 
circumstellar matter 



1 INTRODUCTION 

Symbiotic stars are generally believed to be binary systems 
of a mass losing giant and a hot star, which is usually a 
white dwarf (Kenyon 1986). The giant component loses a 
large amount of mass in the form of a slow stellar wind, a 
significant fraction of which may be ionized by the strong 
UV radiation emanating from the hot component (e.g. Tay- 
lor & Seaquist 1984). The physical extent of the ionized 
region should be sensitively dependent on the strength of 
the ionizing radiation and the mass loss rate of the giant 
component. In this respect, studies of ionization structures 
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of symbiotic stars are very important in understanding the 
mass loss process in symbiotic stars (e.g. Skopal 2003). 

Symbiotic stars exhibit very unique spectroscopic fea- 
tures, which are Raman scattered lines by neutral hydrogen. 
The emission features around 6825 A and 7082 A appearing 
in more than half of the symbiotic stars are Raman scat- 
tered O VI 1032, 1038 doublets first identified by Schmid 
(1989). The Raman scattering origin for these features is 
strongly supported by many spectropolarimetric observa- 
tions (e.g. Schmid & Schild 1994, Schmid et al. 2000, Harries 
& Howarth 1996) and contemporaneous far UV observations 
(e.g. Espey et al. 1995, Birriel, Espey, Schulte-Ladbeck 1998, 
2000). The Raman scattering cross sections for O VI AA 
1032,1038 are of order 10 -22 cm 2 . Therefore, the existence 
of these features implies that a non- negligible fraction of the 
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scattering region is characterized by a very high H I column 
densities Nhi ^ 10 22 cm~ 2 . 

Van Groningen (1993) found other Raman scattered 
features including Raman scattered He II lines in the sym- 
biotic star RR Telescopii. In particular, Raman scattered 
He II line at around 4850 A was found in the spectrum of 
the young planetary nebula NGC 7027 by Pequignot(1997). 
Because a He II ion has the same atomic structure with 
four times large level spacings as a hydrogen atom does, the 
emission lines associated with the energy levels with even 
principal quantum numbers of He II have a slightly shorter 
wavelengths than those of H I emission lines. This implies 
that the Raman scattering cross sections for these He II 
lines are so high that the scattering processes can operate in 
broader range of objects including young planetary nebulae 
(Nussbaumer, Schmid & Vogel 1989). 

In the same line of reasoning, UV continuum radiation 
around Ly/3 will be Raman scattered to be redistributed 
around Ha, forming broad wings. In this case, the wing 
profiles reflect the wavelength dependence of the scatter- 
ing cross section, which is approximately inversely propor- 
tional to the squared wavelength difference of the incident 
radiation and the Ly/3 line center. These broad Ha wings 
with profiles well-fitted by a profile oc AA~ 2 were found 
in objects including many symbiotic stars (Lee 2000), some 
young planetary nebulae such as IC 4997 (Lee & Hyung 
2000) and M2-9 (Balick 1989), and post AGB stars (van 
de Steene, Wood, & van Hoof 2000). Recently Arrieta & 
Torres-Peimbert (2003) added 12 more objects that are also 
in late stages of stellar evolution. Considerable mass loss is a 
common feature for all of these objects, where one can find 
a large amount of neutral matter surrounding the central 
hot stars. It is also interesting that a large fraction of these 
objects exhibit bipolar nebular morphology (e.g. Schwarz & 
Monteiro 2003, Schwarz & Corradi 1992). 

Lee et al. (2003) used the 3.6 m Canada-France-Hawaii 
telescope to perform spectroscopy of the symbiotic star 
V1016 Cygni. They found the Raman scattered He II 6545 
that is converted from He II A 1025 emission line. They 
adopted a simple scattering geometry to compute the broad 
Ha wing template profiles and isolate the 6545 feature, from 
which they were able to put a lower limit of the mass loss 
rate of the Mira variable. However, the strengths of Raman 
scattered features are determined by a combination of the 
covering factor and neutral hydrogen column density (Nhi) 
of the scattering region for given incident UV radiation. The 
degeneracy of the two factors should be lifted in order to put 
more stringent constraints on the mass loss processes occur- 
ring in symbiotic stars. 

There exists conspicuous asymmetry in the branching 
ratio of Raman scattering relative to Rayleigh scattering 
around Lya, where the branching ratio is an increasing func- 
tion of the incident wavelength (Yoo, Bak, & Lee 2002). 
Therefore, continuum radiation redward of Ly/3 will be more 
easily Raman scattered than continuum blueward of Ly/3. 
This tendency increases for higher scattering optical depths, 
which may lead to a redward shift of the center location 
of Ha wings as Nhi increases. Furthermore, the Raman 
conversion rate increases in the wavelength region, where 
He II A 1025 is located. Therefore, if He II A 1025 exhibits a 
symmetric emission profile with the peak at its line center, 
the peak position of the Raman scattered He II 6545 will 



move toward Ha, where the amount of shift will be mainly 
determined by the neutral hydrogen column density. 

In this paper, we perform Monte Carlo simulations to 
compute the center shifts of the Ha wings and the Raman 
scattered He II 6545 features for various Nhi, which may 
be useful in lifting the degeneracy of the column density 
and the covering factor that determine the strengths of the 
Raman scattered features. In the following section we briefly 
describe the basic ideas of the center shifts. The numerical 
results from our Monte Carlo calculations are presented in 
Section 3. Finally in Section 4 we discuss our results and 
observational ramifications. 



2 MODELS AND MONTE CARLO 
PROCEDURES 

In a symbiotic star, the gravity of the white dwarf signif- 
icantly affects the stellar wind flow from the giant, as is 
shown by the SPH computations by Mastrodemos & Morris 
(1998). However, it will still hold that the matter density 
will decrease as the distance from the giant center increases. 
As a first approximation, the ionization structure was in- 
vestigated by Taylor & Seaquist (1984), assuming that the 
circumstellar distribution is radial. With this assumption, 
the ionized region will take a conical shape, which is deter- 
mined by the relative strength of H ionizing UV radiation 
Lh and the mass loss rate M of the giant. They introduced 
a parameter X H + defined as 

Xh+ = a B (H,T e )a H DLH l^J W 

Here, cib(H, T e ) is the recombination coefficient of hy- 
drogen in a region with electron temperature T e , D is the 
binary separation, /i is the mean molecular weight and an 
is the relative hydrogen particle density. The Stromgren 
boundary is formed where the parameter X H + exceeds the 
value of the function f(6) defined by 

where the angle 6 is the opening angle of the ionized conical 
region (Seaquist, Taylor & Button 1984). 

The Raman scattering optical depth of unity for UV 
radiation corresponding to this spectral region is obtained 
when the neutral hydrogen column density exceeds Nhi ~ 
10 20 cm" 2 (e.g. Nussbaumer et al. 1989, Lee 2000). There- 
fore, the Raman scattered He II 6545 and nearby Ha wings 
are formed in the neutral region, because the neutral column 
density of the ionized region does not exceed the column den- 
sity Nhi = 10 18 cm -2 , for which the Lyman limit radiation 
becomes optically thick. The neutral column density Nhi 
for a given line of sight starting from the UV emission re- 
gion through the neutral scattering region will be inversely 
proportional to the impact parameter with respect to the gi- 
ant. Therefore, the determination of Nhi may put a strong 
constraint on the mass loss rate of the giant component in 
a symbiotic star. 

In this paper, we are mainly concerned with the Raman 
scattered He II 6545 feature and broad Ha wings found in 
symbiotic stars. We adopt a view that the broad Ha wings 
are formed through Raman scattering of continuum around 
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Figure 1. A spectrum around He II A 1085 of the symbiotic star 
RR Telescopii obtained with the ORFEUS. The He II A 1085 
line is well fitted by a single Gaussian fx oc /oe~( A ~ A °) / AA , 
where AA = 0.1 A and f = 4.0 X 10~ 12 erg s _1 cm" 2 A -1 . The 
continuum is shown by a solid line and dotted lines represent 2.5<r, 
where the statistics were obtained from the wavelength interval. 



Figure 2. The Raman scattering cross section overplotted to 
the He II 1025 emission line given by a single Gaussian with 
AA = 0.08 A. Note that the cross section increases steeply toward 
the Ly/3 line center, which will affect the center location of the 
Raman scattered He II 6545 feature. 



Ly/3, as is proposed by Nussbaumer et al. (1989). The rel- 
ative strength of Raman scattered He II 6545 feature with 
respect to that of Ha wings are determined basically by the 
equivalent width of He IIA 1025 and the scattering optical 
depth. 

Schmid et al. (1999) studied O VI AA 1032, 1038 lines 
using ORFEUS (Orbiting and Retrievable Far and Extreme 
Ultraviolet Spectrometers). In Fig. 1, we show a spectrum 
around He II A 1085 of the symbiotic star RR Telescopii 
obtained with the ORFEUS, in order to obtain a hint of 
He II A 1025 which suffers heavy interstellar extinction 
due to closeness to H I Ly/3. It is well fitted by a single 
Gaussian fx oc / e" (A " Ao)2/AA2 , where AA = 0.1 A and 
fo = 4.0 x 10~ 12 erg s _1 cm~ 2 A -1 . In order to estimate 
the equivalent width of He II A 1085, we fit the continuum 
and the He II A 1025 emission line by a solid line in the fig- 
ure, from which the continuum flux density is chosen to be 
f c — 4 x 10~ 13 erg cm -2 s _1 A -1 . By dotted lines we show 
the continuum levels higher and lower by 2.5<r. We obtain 
the equivalent width 3.7 A, and put a lower bound 2.0 A 
from this measurement. Using the Case B recombination 
theory we may infer the equivalent width of He A 1025 with 
the additional assumption that the continuum level does not 
change much. In the current work, we adopt a value 0.5 A, 



which appears to fit to the spectrum of V1016 Cyg obtained 
by Lee et al. (2003). 

When the scattering optical depth is small, the Ra- 
man conversion efficiency will be approximately propor- 
tional to the cross section that is again approximately given 
by oc (A — Aiy/3) -2 (e.g. Lee & Lee 1997, Lee & Hyung 
2000). Therefore, the peak position of the Raman feature is 
determined by the location of the maximum of the function 



P(A) 



p e 



-l(X-X Ha+ )/AX] 2 



(3) 



where \ He + and \hyi3 are the center wavelengths of 
He II A 1025 and Ly/3, respectively. Here, AA = 0.08 A is 
the width of the single Gaussian for incident He II A 1025 
emission line, which was chosen in this work. Here, we de- 
fine Y = \h y j) — ^He+j an d m e and m p are masses of an 
electron and a proton, respectively. In this work, we choose 
X Lyl3 = 1025.7222 A and \ He + = 1025.273 A. In Fig. 2, we 
show the Raman scattering cross section and the He II 1025 
emission line given by a single Gaussian with AA = 0.08 A. 

By performing a simple differentiation, we find 
Y - VY 2 - 4AA 2 



(4) 
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where \ m ax,i is the wavelength around X He +, at which the 
Raman scattered He II will exhibit the peak. In the limit 
A A <C Y , the peak of Raman conversion rate is found at 

\ ma *,i ^ A He+ + AA 2 /y ~ \ He+ + 0.016 A, (5) 

for our assumed value AA = 0.08 A. It is interesting to note 
that the amount of peak shift is proportional to the squared 
width of the incident radiation with a Gaussian profile. This 
shift is translated into the redward shift of Raman scattered 
He II 6545 feature by an amount 

\ max ,o m A ffe + flam +0.016x6.4 2 A = A ffe+iflom +0.66 A, (6) 

which is quite significant and should be easily observable. 

As the scattering optical depth increases, a significant 
fraction of Raman scattered photons are formed after a few 
Rayleigh scatterings, resulting in a constant Raman conver- 
sion rate, under the assumption of the absence of absorbing 
opacities such as dust and little opacity near Ha (e.g. Schmid 
1996, Lee & Lee 1997). In this highly optically thick limit, no 
shift in peak position will be seen. This implies that the red- 
ward shift of Raman scattered feature He II 6545 decreases 
as the neutral hydrogen column density increases, render- 
ing the amount of redward shift a useful indicator of Nhi- 
The multiple scattering effects are difficult to handle ana- 
lytically, and we adopt a Monte Carlo approach to obtain 
the quantitative relation between the peak shift and Nhi- 

Monte Carlo techniques have been applied usefully to 
understand basic properties of Rayleigh and Raman scat- 
tered radiation and interpret various spectroscopic and po- 
larimetric observational data (e.g. Schmid 1995, 1996). In 
this work, we use the Monte Carlo code developed by Yoo, 
Bak & Lee (2002), in which they adopted the branching ratio 
rt of Raman and Rayleigh scattering cross sections around 
Ly/3 given by 

tt = 0.1342 + 12.50 ( A ; XLyP ) . (7) 

This situation is illustrated in Fig. 1 of Yoo et al. (2002), 
where the branching ratio was represented by a dotted line. 
Although it is highly unrealistic to characterize the scatter- 
ing region by a single value of Nhi, we adopt the scattering 
regions in this way in order to find the amount of center 
shift dependent on Nhi- 

The Raman scattering branching ratio increases as the 
wavelength increases fairly steeply according to Yoo et al. 
(2002). For higher column densities the Raman wings will 
extend to broader range of wavelengths, where the difference 
in the branching ratios redward and blueward of the line 
center will increase significantly, strengthening the red wing 
part selectively. This effect may lead to the redward shift of 
the Ha wing center as Nhi increases. 



3 RESULTS 

3.1 Center Shift of Ha Wings 

In Fig. 3, we show Ha wings produced by our Monte Carlo 
code for Raman scattering of continuum radiation around 
Ly/3 for various neutral hydrogen column densities in order 
to determine the wing center. The vertical axis represents 
the Raman conversion rate, which is the ratio of the number 
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Figure 3. Ha wings produced by our Monte Carlo code for 
Raman scattering of continuum radiation around Ly/3 for various 
neutral hydrogen column densities in order to determine the wing 
center. The vertical axis represents the Raman conversion rate f\. 
The dashed line shows the best fitting profile with the functional 
dependence oc AA~ 2 = (A — Aq)~ 2 in the range fx < 0.5. 



of incident UV photons to that converted into optical ones. 
As is explained by Lee & Hyung (2000), when the medium 
becomes optically thick, the conversion rate saturates to a 
constant value near 0.6. 

A least square method was adopted to locate the wing 
center. Using the fact that the far wing parts are fairly 
well approximated by profiles oc AA~ 2 = (A — A () )~ 2 , by 
varying Ao we obtained a AA -2 profile that best fits each 
Monte Carlo result for those parts with the conversion 
rate less than 0.5. For illustrative purpose, one such fit for 
Nhi = 10 21 cm~ 2 is shown in Fig. 3 by a long-dashed line. 
Even though (A — \hu)~' 2 profiles provide overall satisfac- 
tory fits to simulated Ha profiles, it should be noted that 
small deviations are also apparent especially near Ha center. 

In Fig. 4, we plot the center locations as a function 
of Nhi- The result shows the redward wing shift of Ha 
wings as Nhi increases. However, the center shift amounts 
to only 20 km s _1 for a change of Nhi from 10 20 cm~ 2 to 
10 21 cm~ 2 . Furthermore, the exact center location is am- 
biguous as is shown by the considerable scatter in Fig. 4, 
which is mainly due to the statistical nature of our Monte 
Carlo calculations. From this simulation, we may conclude 
that the center shift of Ha wings in this range of Nhi does 
not have much observational significance. 
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Figure 4. Center locations of Ha wings as a function of Nhi- 
The scatter shows the statistical nature of the Monte Carlo cal- 
culations. The center shift amounts only 20 km s" 1 for a change 
of N HI from 10 20 cm~ 2 to 10 21 cm" 2 . 



3.2 Raman Scattered He II 6545 

Assuming a similar nebula condition discussed in Section 2, 
we prepare He II A 1025 line as an input spectrum repre- 
sented by a single Gaussian with a width of AA = 0.08 A 
shown in Fig. 2. As is discussed later, this choice explains 
fairly well the spectrum around Ha of the symbiotic star 
V1016 Cyg obtained with the 3.6 m Canada-France-Hawaii 
Telescope in 2002 (e.g. Lee et al. 2003). First, we investi- 
gate the Raman scattered He II 6545 feature without any 
contribution from the UV continuum. 

In Fig. 5, we show the profiles of the Raman scattered 
He II 6545 from a scattering region with neutral hydrogen 
column densities Nhi = 2 x 10 19 , 8 x 10 19 , 4 x 10 20 , and 9.5 x 
10 20 cm" 2 . In order to clarify the shift of the peak position, 
by long dashed lines we show the single Gaussian line profile 
of the Raman scattered He II 6545, which would appear in 
the case where the Raman conversion rate is constant. The 
Monte Carlo results are shown by solid lines. In this figure, 
we immediately see that the center location moves redward 
or toward Ha. The result for Nhi — 2 x 10 19 cm" 2 is quite 
noisy due to low photon statistics. 

A least square method was applied to the Monte Carlo 
results in order to obtain the fits and locate the peak po- 
sitions of the Raman scattered features. We note that all 
the Raman features are fitted satisfactorily by a Gaussian 
function with a width AA 6545 = 6.4 2 x 0.08 A = 3.3 A. 

In Fig. 6, we plot the amount of center shift as a func- 
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Figure 5. Profiles of the Raman scattered He II 6545 from a 
scattering region with N HI = 2 X 10 19 ,8 X 10 19 ,4 X 10 20 , and 
9.5 X 10 20 cm" 2 . The incident He II 1025 emission line is assumed 
to be given by a single Gaussian shown in Fig. 2. Long dashed 
lines show Gaussian line profiles of He II 6545 that would be 
obtained for a constant Raman conversion rate. Note that the 
peak locations move redward by different amounts for various 
Nhi- 



tion of Nhi- For a column density Nhi 



10 2 



the 



peak position shift amounts to AAo ~ 0.9 A, which is quite 
significant and similar to the value discussed in the previ- 
ous section. The amount of center shift decreases as Nhi 
increases in the region Nhi > 5 x 10 19 cm" 2 , and no center 
shift is seen for Nhi > 10 22 cm" 2 . At this high column den- 
sity, the Raman conversion rate for He II A 1025 saturates to 
a constant value ~ 0.6, which explains no shift of the peak 
position. The decrease of the peak shift as Nhi indicates 
the decrease of the slope of the Raman conversion rate as a 
function of the scattering optical depth, which approaches 
zero in the optically thick limit. 

For lower column densities, the peak shifts are quite 
large exceeding 1 A. This implies that in an optically thin 
limit, the Raman conversion rate is not given by a simple 
relation that is proportional to (X — Xl v 0)~ 2 ■ It appears that 
more redward peak shift in low column densities may be 
explained by introducing a parameter a so that the Raman 
conversion rate is described in terms of this parameter by 



Pi(A) 



Poe 



-l(X-X Hc+ )/AX] 2 



(A — Xlvp) 2 



(8) 



where the peak location coincides with one of the roots of 
the derivative p'i(A). The introduction of a new parameter a 




Figure 6. Center shift of the Raman scattered He II 6545 feature 
as a function of Nhi- The amount of center shift decreases as Nhi 
increases, and no center shift is seen for Njji > 10 22 cm~ 2 . 

may be justified by the following argument. In an optically 
thin scattering region characterized by the total scattering 
optical depth r and the Raman branching ratio r b , the Ra- 
man conversion efficiency including multiple scattering ef- 
fects may be approximated by an infinite geometric series 

C Ram ~ rr b + t(1 - r b )rr b + [r(l - r b )] 2 rr b H 



Since we approximately have r oc (A — Xlvp)' 2 , this leads to 
a result in the functional form Cn a m oc [(A — XLyp) 2 — a] -1 . 
With the assumption of an optically thin medium, it follows 
that a < Y 2 . 

By differentiating pi(A) and defining 8y = A — Xl v i3, we 
obtain a cubic equation in 5y, which is 

(8y) :i + Y{Syf + (AA 2 - a)(Sy) - aY = 0. (10) 

Upon investigating the roots of the cubic equation for vari- 
ous values of a, we find that when a ~ 0.09 A , the peak shift 
can be as large as 1.2 A found for Nhi = 0.95 x 10 19 cm~ 2 
in our Monte Carlo calculation. 

In the cases of optically thick scattering regions, the 
Raman conversion rate depends on Nhi in a non-linear way, 
which gives rise to the various amount of the peak shift for 
Nhi- In turn, the nonlinearity is attributed to the multiple 
Rayleigh scattering effect before Raman conversion takes 
place. 

In Fig. 7, we show combined spectra of UV continuum 



Figure 7. Combined spectra of UV continuum and He II A 1025 
that are Raman scattered from scattering regions with various 
Nhi- The strengths of both the H« wings and the He II A 6545 
feature vary for different values of Nhi, or the scattering optical 
depth. Due to the unknown covering factor of the scattering re- 
gion, these simulated profiles should be multiplied by a constant 
factor to fit an observed spectrum. 

and He II A 1025 for various Nhi ■ The strengths of both the 
Ha wings and the He II A 6545 feature vary as the column 
density of the scattering region mainly due to the difference 
in scattering cross section. However, the observed strength is 
directly proportional to the covering factor of the scattering 
region, and hence these two factors are degenerate. If we fix 
the covering factor and geometrical shape of the scattering 
region, the Raman conversion efficiency is determined by the 
scattering optical depth in a non-linear way. Therefore, with 
the unknown covering factor of the scattering region, these 
simulated profiles should be multiplied by a constant factor 
to fit an observed spectrum. 

In Fig. 8, we show the spectrum of the symbiotic star 
V1016 Cyg obtained with the 3.6 m Canada-France-Hawaii 
Telescope (see Lee et al. 2003 for observation details). We 
overplot the simulated Ha wings and He II A 6545 features 
for various column densities onto this spectrum. It appears 
that the simulated profiles with Nhi in the range consid- 
ered in this figure provide good fits to the observed spectrum 
blueward of the [N II] 6548 line. Therefore, in order to deter- 
mine Nhi, it is necessary to locate the exact peak position, 
which can be done only after a very careful de-blending of 
[N II] A 6548. 

From the fact that He II 6545 feature is observed outside 
the region where the Raman conversion rate is saturated 



Raman Scattered He II in Symbiotic Stars 7 




2000 
1500 
1000 
500 


1.5x10' 
* 10 s 
5x10 s 


4xl0 4 
3xl0 4 
" 8x10" 
10 4 





1 






; . aJ 


\ 

u 


II. 

1 

1— — : 



6550 6600 



6550 6600 

MA) 



6500 



6550 



6600 



Figure 8. A spectrum of the symbiotic star V1016 Cyg obtained 
with the 3.6 m Canada-France-Hawaii Telescope The simulated 
Ha wings and He II 6545 feature are overplotted represented by 
dotted lines. 



and constant, the scattering region should be characterized 
by N HI < 10 21 cm -2 . Noting that Raman scattered O VI 
6825 and 7082 features are formed in a region with Nhi ^ 
10 22 cm -2 , the scattering region responsible for He II 6545 
should include the O VI scattering region and be much more 
extended. 



4 DISCUSSION AND OBSERVATIONAL 
RAMIFICATIONS 

4.1 Electron Scattering Wings 

Electron scattering of Ha line photons is one of physical 
mechanisms that have been proposed to be the origin of the 
broad Ha wings in many celestial objects including active 
galactic nuclei and Wolf-Rayet stars (e.g. Mathis 1970, Wey- 
mann 1970, Hillier 1991). We generated electron scattering 
wings using a similar Monte Carlo technique, which is shown 
in Fig. 9. 

Because the Ha emission line is well-approximated by 
a single Gaussian fx = 1650e" (A_Aff " )/AAff ° , we used this 
Gaussian as an incident radiation. The temperature of the 
electron cloud is chosen to be T e = 2.5 x 10 4 K. The cov- 
ering factor of the electron scattering region is assumed 
to be C e i = 3.5 x 10~ 2 and the electron column density 
iV e = 10 24 cm~ 2 , for which the Thomson optical depth 
r e — 0.67. It is also noted that similar profiles will be ob- 



Figure 9. Electron scattering wings overplotted to the spectrum 
of V1016 Cyg obtained with the CFHT. The covering factor of 
the electron scattering region is assumed to be C e i =0.05 and the 
electron column density N e = 10 24 cm -2 , for which the Thomson 
optical depth r e = 0.67. The temperature of the electron cloud is 
chosen to be T e = 2.5 X 10 4 K. The top panel shows the electron 
scattering wing represented by a solid line, which was fitted us- 
ing the theoretical profile shown by a dotted line investigated by 
Weymann (1970). 



tained by a choice of r e and C e i with their product fixed as 
far as the electron plasma is optically thin. 

In the top panel, we show the electron scattering wing 
produced by a Monte Carlo technique, which is represented 
by a thin solid line. According to Weymann (1970), electron 
scattering wings formed in an optically thin region may be 
approximated by exponential profiles represented by a thick 
solid line with a functional form 



/e(A) = foe 



(A-A ffQ )/AA e 



(11) 

Here, the profile width AA e = 22 A coincides with the 
Doppler width corresponding to the thermal electron veloc- 
ity with T e = 2.5 x 10 4 K, from which the electron tempera- 
ture is strongly constrained. The apparent discrepancy with 
the theoretical profile is attributed to our choice of rather 
high r e and the fact that the incident Ha emission profile 
is not a Dirac delta function but a Gaussian with a finite 
width. 

In the middle and bottom panels, the Monte Carlo re- 
sults are represented by dotted lines and the solid lines show 
the CFHT spectrum of V1016 Cyg. The fit is quite satis- 
factory. However, the Raman scattering origin of Ha wings 
naturally explains the existence of the He II 1025 feature. 
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A more stringent constraint may be put from careful analy- 
ses of the broad H/3 wings and He II 4850 feature, because 
Raman scattering cross sections and branching ratios dif- 
fer from those corresponding to features around Ha. Future 
spectroscopic observations covering both Ha and H/3 regions 
simultaneously may shed much light on this point. 

4.2 Observational Ramifications 

An interesting possibility is that if the He II emission is not 
isotropic but the continuum is, then the equivalent width 
of He II A 1025 measured by a hypothetical observer lo- 
cated in the neutral scattering region does not have to be 
equal to that measured by an observer on Earth. Accord- 
ing to Mastrodemos & Morris (1998), some fraction of the 
stellar wind from the giant component may be captured by 
the white dwarf component in symbiotic stars to form an 
accretion disk. If we accept this picture, the line emission 
region may take a disk-like shape, where more line radiation 
is expected to the normal direction of the orbital plane or 
the accretion disk. Using the ORFEUS data, Schmid et al. 
(1999) pointed out that O VI AA 1032, 1038 doublets exhibit 
different profiles from those of their Raman scattered coun- 
terparts, which appears to support this picture. In this case, 
for an observer with a large inclination to the orbital plane, 
the equivalent width that is estimated from the Raman scat- 
tered He II 6545 will be smaller than that estimated from 
the far UV spectroscopy performed for the same object. 

In the work of Lee et al. (2003), the neutral column 
density Nhi was assumed to be the value that gives a unit 
scattering optical depth, of which the validity has not been 
pursued in a rigorous way. Therefore, if this choice of Nhi 
is lower than the true value, then the covering factor of the 
scattering region should be larger than they estimated, and 
vice versa. However, in order to locate the peak position very 
accurately, we have to subtract [N II] A 6548 feature reliably 
using the 3 times stronger [N II] A 6584, which requires very 
careful data reductions. [N II] lines are weak or absent in the 
higher nebular density S-type systems, for which case the 
line center of the Raman scattered He IIA 6545 can poten- 
tially be determined without this complication. As is shown 
by this study, an accurate value of the mass loss rate of the 
giant can be obtained by a more realistic model equipped 
with a range of Nhi with independent information about 
the equivalent width of He II A 1025. 

It is particularly notable that the continuum photon 
number flux is diluted by a factor of (Ai/A ) 2 due to the 
incoherency of the Raman scattering process. Because the 
equivalent width is a measure of the line photon number 
flux normalized by the continuum photon number flux, the 
computation of the equivalent width of the He II A 1025 
presented by Lee et al. (2003) is erroneous and overestimated 
by a factor of 6.4. The correct equivalent width should be 
about 0.51 A which is similar to the value adopted in the 
current work. 

It is very interesting to note that the asymmetric 
branching ratio of the Raman scattering cross section gives 
rise to asymmetric wing profiles with the stronger red part 
around Ha. Regarding this, spectropolarimetric studies of 
Raman scattered O VI 6825 and 7082 features by Harries 
& Howarth (1996) show the flip of polarization direction 
occurring in the far red part, implying the existence of a 



third emission (or scattering) component moving away from 
the main scattering region (or emission region). Therefore, 
when examining this kind of component that is moving away, 
the effect from the asymmetry of atomic physics should be 
corrected first. More spectropolarimetry will be very useful 
to investigate the existence of similar phenomenon in Ha 
wings. 

Less scattering cross sections and more branching possi- 
bilities may lead to more diverse features in Raman scattered 
lines around H/3 and higher Balmer series lines. Therefore, 
simultaneous spectroscopy for the Raman scattered features 
around H/3 and higher Balmer series lines will shed more 
light on the mass loss processes and the ionization struc- 
tures of symbiotic stars. Similar center shift phenomena are 
expected to Raman scattered He II features formed blue- 
ward of higher Balmer series lines. Because center shifts are 
also sensitive to relative kinematics between the UV emis- 
sion source region and the scattering region, careful data 
analyses will be required. We conclude that Raman features 
provide very important and interesting diagnostic tools for 
investigating the mass loss processes and ionization struc- 
tures in symbiotic stars. 
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